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Procede de localisation d'un element d'interet contenu dans un 
objet tridimensionnel, en particulier lors d'un examen de 
stereotaxic en mammographie. . : ^ ; 

L'invention conceme la localisation d'un element d'interet contenu 
dans un objet tridimensionnel h partir des positions de regions 
d'interet homologues correspondant audit 616rnent d'interSt et 
apparaissant dans un ensemble damages stereotaxiques dudit objet. 

Bien que l'invention s'applique k tout ensemble d'images 
stereotaxiques d'un objet tridimensionnel quelconque, elle est 
particulierement utile dans le domaine medical et notamment lors dun 
examen de stereotaxic en mammographie pour la localisation de 
microcalcifications dans un sein. 

Un examen de stereotaxic en general permet d'acceder avec une 
grande precision, de l'ordre du millimetre, a un point donne d'un objet 
tridimensionnel, k partir de deux projections bidimensionnelles de cet 
objet acquises par exemple selon deux incidences angulaires opposees 
de part et d'autre de la normale au plan du recepteur damages delivrant 
lesdi tes proj ections . 

A partir de la connaissance, d'une part, de la position de la 
projection de ce point dans le plan de chaque image stereotaxique 
obtenue, et, d'autre part, de la geometrie de Tappareil de prises de 
vues stereotaxiques ay ant conduit a la construction des deux images, il 
est possible de calculer la position exacte spatiale de ce point dans 
l'objet tridimensionnel, par calcul trigonometrique. 

En resume, un examen stereotaxique necessite : 

de faire au moins deux images de l'objet a differentes angulations, 
de connaitre parfaitement la geometrie d'acquisition du systeme 
stereotaxique, et 

de pouvoir localiser sur les differentes images stereotaxiques 
obtenues, le lieu de la projection de Tel^ment d'interet choisi. 

Les deux premiers points ne posent pas de probl&me majeur. 

Actuellement, en mammographie, 1'examen de sterdotaxie est 
effectue & l'aide d'un mammographe equip6 d'un dispositif de prises de 
vues stereotaxiques. Le mammographe comprend un tube & rayons X, 



situe k rextremit6 d'un premier bras qui est mobile autour d'un axe et 
emettant un rayonnement X vers un recepteur situ6 k rextremite d f un 
autre bras. Entre le tube et le recepteur sont disposes, d'une part, un 
plateau porte-sein, ou support patient, et, d'autre part, un plateau de 
compression assurant le rnaintien en place du sein lors des 
mammographies. Le recepteur d'images peut Stre un recepteur 
numerique, tel qu'une camera CCD par exemple, d61ivrant des images 
stereotaxiques numerisees. Les prises de vues st6r6otaxiques 
necessitent la rotation du tube a rayons X autour des plateaux porte- 
sein et de compression, selon deux orientations successives opposees 
de part et d'autre de sa position initiale perpendiculaire au plan du 
recepteur damages. 

A partir des images stereotaxiques numerisees obtenues, la 
localisation de l'element d'interet sur les differentes images est une 
operation delicate, en particulier en marnmographie, puisqu'elle fait 
appel k une mise en correspondance, ou appariement ("matching" en 
langue anglaise) de regions d'interet homologues correspondant, dans 
deux images stereotaxiques differentes, k un meme element d'interet 
situe dans le sein. 

Ainsi, en marnmographie, on cherche a apparier sur les differentes 
images stereotaxiques obtenues, les regions d'interet qui correspondent 
a une meme microcalcification. 

Or, l'appariement des microcalcifications projetees est un 
probleme difficile car les microcalcifications ne se ressemblent pas 
d'une image a l'autre. Leur forme et leur contraste peuvent changer 
ainsi que leur agencement dans l'espace. Elles peuvent se superposer 
avec des zones fibreuses ou avec d'autres microcalcifications. 

Jusqu'a present, Tappariement etait effectue entierement 
manuellement par le radiologue. Ce dernier choisissait quelles etaient 
les regions d'interet en correspondance sur les deux images 
stereotaxiques angulees, puis rentrait les coordonnees 
bidimensionnelles de ces regions dans un calculateur qui donnait la 
position spatiale de la microcalcification, c'est-k-dire les coordonnees 
tridimensionnelles de cette demiere. A partir de ces coordonnees, le 
radiologue pouvait ajuster par exemple la position d'un porte-aiguille 



de fa§on k ce que Taiguille vienne k Taplomb de la microcalcification 
pour permettre sa ponction par exemple. 

Or, cet appariemment entitlement manuel s'av&re etre une etape 
longue et parfois imprecise. 

L'invention vise k apporter une solution k ce problfeme. 

L'invention a pour but de reduire sensiblement la dur6e d'un 
examen st6r6otaxique et d'obtenir une meilleure precision 
d'appariement. 

D'une fagon tres generate, selon l'invention, on definit une mesure 
de similarite permettant d'evaluer la ressemblance entre une region 
d'interet (microcalcification par exemple) selectionnee (manuellement 
par le radiologue par exemple, ou bien automatiquement par un 
procede de selection automatique) sur une des images stereotaxiques 
angulees et une region d'interet d'une deuxieme image stereotaxique. 

A partir de cette mesure de similarite, on selectionne, pour obtenir 
la region d'interet homologue, un certain nombre de "candidats" 
appartenant a une zone de recherche et ayant une forte ressemblance 
avec la microcalcification de depart. 

Ensuite, en utilisant les principes de la stereotaxic, on calcule les 
coordonnees tridimensionnelles de ces candidats, puis on les reprojette 
sur Timage de centrage. On calcule k nouveau la similarite entre ces 
points reprojetes et la microcalcification initialement choisie. En 
effet, si le candidat correspond effectivement k la microcalcification 
de depart, alors sa projection calculee correspondra effectivement k 
Timage de cette microcalcification sur l'image de centrage. 

Enfin, a Taide de ces deux valeurs de similarite calculees pour 
chaque candidat, on propose le meilleur candidat. 

En theorie, on pourrait s'abstenir de Tetape de reprojection 
utilisant Timage de centrage et se contenter de Tappariement tel qu'il 
vient d'etre expose pour selectionner la region d'interet homologue de 
celle selectionnee sur Timage cible. Cependant, il est nettement 
preferable d'utiliser la reprojection des candidats sur Timage de 
centrage pour verifier la premi&re s61ection effectu6e, notamment dans 
le domaine de la mammographie ou les microcalcifications peuvent 
changer de forme et de contraste d'une image k Tautre. 



En d'autres termes, l'invention propose tout d'abord un proced6 de 
localisation d'un element dlnterdt contenu dans un objet 
tridimensional k partir des positions de regions d'interets homologues 
correspondant audit element d'interet et apparaissant dans un ensemble 
d'images stereotaxiques dudit objet. Ce precede comporte une etape de 
selection dans une premiere image stereotaxique, ou image cible, d'une 
premiere region d'interet, en particulier une microcalcification, dite 
region d'interet cible. Le precede comporte egalement un appariement 
de ladite premiere region avec une deuxieme region d'interet 
homologue de la premiere et apparaissant dans une deuxieme image 
stereotaxique. 

Selon une caracteristique generate de l'invention, les images 
stereotaxiques etant numerisees, on selectionne dans Tetape de 
selection un pixel cible dans ladite region d'interet cible. Dans Tetape 
d'appariement, on genere autour du pixel cible selectionne, une fenetre 
cible de caracteristiques dimensionnelles choisies et contenant ladite 
region d'interet cible. On determine alors dans la deuxi&me image, 
selon un critere de selection predetermine, un ensemble de pixels et on 
genere autour de chaque pixel selectionne une deuxieme fenetre de 
memes caracteristiques dimensionnelles que ladite fenetre cible. On 
effectue un traitement de correlation entre les niveaux de gris des 
pixels de chaque deuxieme fenetre et les niveaux de gris des pixels de 
la fenetre cible de fa?on a obtenir une valeur de correlation pour 
chaque deuxieme fenetre. On identifie alors la region dlnteret 
homologue de la region d'interet cible h partir de Tanalyse de 
Tensemble des valeurs de correlation ainsi obtenues. On minimise 
ainsi les risques d'erreur d'appariement entre des regions dlnteret 
homologues. 

Selon un mode de mise en oeuvre, Tanalyse des valeurs de 
correlation obtenues comporte la selection d'un nombre choisi de 
maxima ou de minima de correlation, la selection de la region 
homologue s'effectuant parmi celles dont la valeur de correlation 
associee est Tun de ces maxima ou minima de correlation. 

En fait suivant que Timage est "normale" ou inversee (une image 
"normale" etant definie comme presentant des niveaux de gris et des 



contrastes qui sont semblables & ceux habituellement visibles sur des 
cliches de film) on utilisera respectivement les maxima ou les minima. 

On peut egalement prevoir une determination de la dynamique des 
maxima ou minima selectionnes, et la comparaison des valeurs de 
dynamique obtenues k un seuil. , 

Le traitement de correlation utilisd peut comporter un traitement 
de correlation normee, centree ou non, ou bien un traitement de 
difference normee. 

En variante a Tappariement evoqu6 ci-dessus, Tinvention propose 
egalement un appariement base sur Telaboration de caracteristiques 
predeterminees de toute region d'interet, en particulier des 
caracteristiques physiques comme l'elongation, le contraste, le 
gradient, puis une comparaison entre ces listes de caracteristiques 
associees h differentes regions d'interSts. 

Plus precisement, selon une caracteristique generate de 
Tinvention, on selectionne dans T6tape de selection, un pixel cible 
dans ladite region d'interet cible et, dans Tetape d'appariement, on 
genere, autour du pixel cible selectionne, une fenetre cible de 
caracteristiques dimensionnelles choisies et contenant ladite region 
d'interet cible. On determine pour ladite region d'interet cible, k partir 
des valeurs de niveau de gris des pixels de ladite fenetre cible, un 
premier ensemble de valeurs numeriques respectivement 
representatives de caracteristiques predeterminees caracterisant toute 
region d'interet (par exemple des caracteristiques de forme de 
contraste, de gradient ... ) et on stocke ce premier ensemble de valeurs 
numeriques. On determine dans la deuxieme image, selon un critere de 
selection predetermine, un ensemble de pixels appartenant 
respectivement a des regions d'interets dites candidates apparaissant 
dans la deuxieme image et on genere autour de chacun de ces pixels 
selectionnes une deuxieme fenetre contenant ladite region d'interet 
candidate correspondante. On d6termine, pour chaque region d'interet 
candidate h partir des valeurs de niveau de gris des pixels de la 
deuxieme fenetre associee, un deuxieme ensemble de valeurs 
numeriques respectivement representatives desdites caracteristiques 
predeterminees et on stocke ce deuxieme ensemble de valeurs 
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numeriques. On identifie alors la region d'interdt homologue de la 
region d'interSt cible k partir d'un traitement de cornparaison entre le 
premier ensemble de valeurs numeriques et chaque deuxieme ensemble 
de valeurs numeriques. On minimise egalement ainsi les risques 
5 d'erreur d'appariement entre des regions homologues. 

Le traitement de cornparaison peut comporter un traitement de 
minimisation de distance entre les deux ensembles de valeurs 
numeriques. 

Dans les deux variantes d'appariement qui viennent d'etre 
10 evoquees, et.afin de minimiser encore le temps de calcul, on determine 

dans la deuxieme image stereotaxique, une zone epipolaire contenant 
au moins le segment epipolaire relatif au pixel cible. Les pixels de la 
zone epipolaire contiennent alors les pixels selectionnes. Pour reduire 
encore le temps de calcul, au lieu d'effectuer 1'appariement sur tous les 
15 pixels de la zone epipolaire determinee, on selectionne 

avantageusement les pixels de la deuxieme image parmi ceux 
presentant des maxima ou des minima de niveau de gris dont la 
dynamique est superieure a un seuil predetermine. 

En outre, pour s'affranchir des problemes de fond d'image, et pour 
20 minimiser encore les risques d'erreurs d'appariement entre deux 

regions d'interets homologues, il est particulierement avantageux que 
l'appariement selon l'invention comporte prealablement un filtrage des 
images stereotaxiques, par exemple du type dit "chapeau haut de 
forme". 

25 L'utilisation de l'etape de reprojection des elements d'interets 

candidats dans Timage de centrage est independante du type 
d'appariement utilise prealablement. 

Ainsi, selon une caracteristique generale de Tinvention, le procede 
de localisation du type de celui evoque ci-avant cornprend une 

30 selection dans une premiere image stereotaxique num6risee d'une 

premiere region d'interet dite cible, une selection, dans une deuxieme 
image stereotaxique numerisee, k partir d'un premier appariement 
automatique, d'au moins une deuxieme region d'interet susceptible 
d'etre homologue de la region cible, la determination de la position 

35 spatiale d'un pixel candidat d'un element d'interet candidat 



7 



correspondant k ces deux regions d'interSts, la determination de la 
position dans une troisieme image stereotaxique, d'un pixel projete 
correspondant k la projection dans cette troisieme image stereotaxique 
du pixel candidate et un deuxieme appariement automatique entre la 

5 region interet cible et un voisinage du pixel projete definissant une 

region d 'interet projete. 

En pratique, on selectionne avantageusement dans la deuxieme 
image stereotaxique, a partir dudit premier appariement automatique, 
un ensemble de deuxiemes regions susceptibles d'etre homologues de 

10 la region cible et on determine la position spatiale de chaque pixel 

candidat correspondant a chaque paire de regions d'interets form6e par 
la region cible et Tune des deuxiemes regions. On determine 
egalement la position dans la troisieme image de chaque pixel projete 
correspondant et on effectue le deuxi&me appariement automatique 

15 entre la region cible et un voisinage de chaque pixel projete. 

On peut ordonner l'ensemble des deuxi&mes regions obtenues dans 
l'ordre decroissant de leur probability k etre la region homologue de la 
region cible. On peut egalement ordonner l'ensemble des regions 
projetees obtenues dans l'ordre decroissant de leur probabilite k etre la 

20 region homologue de la region cible. On selectionne alors 

avantageusement la region homologue de la region cible comme etant 
celle dont le produit des rangs respectifs dans les deux 
ordonnancements est mininum. 

Le premier appariement automatique peut etre Tun des deux 

25 appariements qui ont ete evoques ci-avant, c'est-^-dire k base d'un 

traitement de correlation ou d'un traitement de comparaison. 

De meme, le deuxieme appariement automatique peut comporter 
un traitement de correlation tel que defini ci-avant, entre le voisinage 
de chaque pixel projete et une fenetre cible contenant la region 

30 d'int6ret cible. Ce deuxieme appariement automatique peut etre 

egalement un traitement de comparaison tel que d6fini ci-avant, entre 
la region d interet cible et chaque region d'interet projet6e. 

D'un point de vue materiel, l'invention s'analyse comme un 
dispositif de localisation d'un element d'interet contenu dans un objet 

35 tridimensionnel k partir des positions de regions d'interets homologues 
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correspondant audit element d'interet et apparaissant dans un ensemble 
d'images stereotaxiques dudit objet. Ce dispositif comporte un 
appareillage de prises de vues stereotaxiques equipe dun support pour 
ledit objet et d'un recepteur d'image tel qu'une camera CCD capable de 
fournir des images stereotaxiques numerisees. II est egalement prevu 
des moyens de visualisation de ces images, par exernple un ecran 
video, tel que celui d'un microordinateur. Des moyens de selection, 
par exernple une "souris" d'un microordinateur permettant de deplacer 
un marqueur sur 1'ecran et de "cliquer" ce marqueur sur un pixel cible 
d'une region d'interet, permettent la selection du pixel cible. II est 
enfin prevu un microprocesseur et des memoires associees pour 
incorporer sous forme logicielle des moyens permettant de realiser 
fonctionnellement toutes les etapes des procedes qui ont ete decrites 
ci-avant. 

D'autres avantages et caracteristiques de l'invention apparaitront a 
l'examen de modes de mise en oeuvre nullement limitatifs, et des 
dessins annexes, sur lesquels : 

la figure 1 est une vue schematique de trois images stereotaxiques 
d'un objet tridimensionnel a partir d'un dispositif de prises de vues 
stereotaxiques; 

la figure 2 illustre tres schematiquement la notion de segment 
epipolaire; 

la figure 3 illustre une zone epipolaire choisie conformement k 
l'invention pour effectuer un appariement; 

la figure 4 est un organigramme general d'un procede selon 
l'invention; 

la figure 5 illustre plus particulierement une etape d'une variante 
du procede utilisant un traitement de correlation; et 

la figure 6 illustre schematiquement plus particulierement un 
procede selon l'invention utilisant une reprojection sur l'image de 
centrage. 

Un examen stereotaxique, en particulier en mammographie, se 
compose d'une serie de trois expositions d'un objet tridimensionnel 1 
(figure 1), par exernple un sein, reposant sur un support 3 et comprime 
par un plateau de compression 2, k l'aide d'un tube k rayons X 
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occupant respectivement trois positions differentes 6, 7 et 8. En 
pratique, on proc&de k une exposition selon un angle de 0° et k deux 
expositions angulees selon deux angles opposes egaux, en pratique 
±15°. 

5 Un recepteur CCD 4 permet d'obtenir ainsi trois images 

stereotaxiques numerisees 10, 11 et 12. Selon une convention 
habituellement utilisee, l'image 10 est l'image de droite, tandis que 
l'image 12 est l'image de gauche et que l'image 11 est l'image de 
centrage. 

10 Un element d'interSt 5, par exemple une microcalcification, 

contenu dans Tobjet tridimensional 1, fournit sur chacune des images 
10, 11 et 12 des regions d'interets respectivement referencees 5d, 5c et 
5g. L'image de centrage sert en particulier k verifier que la 
microcalcification sera accessible durant 1'acte medical eventuel 

15 ulterieur. 

Les coordonnees tridimensionnelles de la microcalcification 5 
peuvent etre obtenues, de fa^on classique et bien connue de l'homme 
du metier, par un calcul simple de trigonometrie a partir de la 
connaissance des coordonnees bidimensionnelles des deux projections 

20 5d et 5g, par exemple. Ces coordonnees sont exprimees par rapport a 

une origine materialisee par le centre d'une croix 9 gravee sur le 
plateau supportant le sein. Et on repere sur chaque image le centre de 
cette croix de reference pour initialiser geometriquement le systeme. 
L'homme du metier sait que, par definition, le segment epipolaire 

25 50 (figure 2) d'un point 5d d'une image 9 d'angulation 01, sur une 

image 12 d'angulation 02, correspond k l'ensemble des points de la 
seconde image 12 pouvant avoir le point 5d comme antecedent. 
Autrement dit, c'est la projection conique par la source d'angulation 
02, du segment de droite 52 joignant ledit point a la source 

30 d'angulation 01. 

Dans le cas ou les images stereotaxiques ont ete obtenues avec des 
angulations telles que les positions du tube k rayons X se situent dans 
un plan parallele aux lignes des images, ce segment epipolaire 50 est 
un segment de droite qui se trouve etre parallele aux lignes de l'image. 

35 Cependant, si Ton n'a pas la connaissance de l'epaisseur de Tobjet 
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tridimensionnel, on ne peut pas determiner la position des extremites 
de ce segment epipolaire. Ce segment epipolaire devient alors une 
droite epipolaire qui, dans la geometrie ^acquisition particuliere 
evoquee ici, se situe sur la meme ligne d'image 12 que le point 
d'interet 5d sur l'image 9. 

En outre, lorsqu'on utilise un detecteur qui presente des 
distorsions g6ometriques, la droite epipolaire se transforme en une 
courbe epipolaire 53 (figure 3). On peut determiner les coordonnees de 
Fensemble des points de la courbe epipolaire sur l'image 12 par une 
calibration geometrique de la camera CCD. II serait alors possible de 
n'effectuer l'appariement selon l'invention que sur l'ensemble de ces 
points. 

Cependant, si pour des raisons de precision de la calibration 
geometrique, on ne peut determiner qu'approximativement les 
coordonnees des points de la courbe epipolaire, on choisira alors 
d'effectuer l'appariement selon l'invention sur un ruban epipolaire 51 
qui encadre la courbe epipolaire. En pratique, on definira une zone 
epipolaire, ou zone de recherche, delimitee par deux lignes paralleles 
passant respectivement par deux points augmentes de part et d'autre du 
point le plus haut et le plus bas de l'epipolaire 53 d'une marge 
predeterminee. Generalement, on opte pour une marge de 3 k 5 pixels. 

On va maintenant decrire plus en detail, en se referant plus 
particulierement aux figures 4 et 5, deux types d'appariement selon 
Tinvention. 

La premiere etape consiste en la selection (etape 400) d'un pixel 
cible 60d au sein d f une region d'interet cible 5d presente par exemple 
dans l'image de droite 10, dite image cible. En pratique, cette 
selection peut s'effectuer manuellement en cliquant un point de la 
region d'interet de Timage concernee apparaissant sur un ecran d'un 
microordinateur. Ceci etant, on peut prevoir egalement un procede de 
selection automatique de pixel. Ce pixel 60d peut etre centre ou non 
dans la region cible. 

On genere ensuite (etape 401) une fenetre cible 6 Id contenant 
ladite region d'interSt cible. Cette fengtre cible peut etre generee 
automatiquement apres avoir applique sur le pixel cible 60d un moyen 
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logiciel classique de croissance de region. En variante, on peut 
egalement prevoir, selon le domaine d'application de l'invention, de 
choisir des fenetres de taille predefinies. Dans tous les cas, on 
obtiendra gen6ralement une fenStre cible rectangulaire ou carree dont 
les c6tes presentent des dimensions de l'ordre d'une dizaine k quelques 
dizaines de pixels. 

Les niveaux de gris de chacun des pixels de cette fenetre cible 
sont ensuite stockes dans une memoire associee k un microprocesseur 
qui incorpore de fa9on logicielle tous les moyens fonctionnels de mise 
en oeuvre du procede decrit. 

A partir de la, on va selectionner dans une etape 402 un ensemble 
de pixels dans une deuxieme image stereotaxique, par exemple 1'image 
de gauche. Cet ensemble de pixels s61ectionnes va constituer une zone 
de recherche 51 k l'interieur de laquelle on va selectionner une ou 
plusieurs regions d'interets candidates susceptibles d'etre homologues 
de la region d'interet cible. 

Comme on l'a explique ci-avant et dans le but de reduire la duree 
d'appariement, cette zone de recherche 51 sera par exemple la zone 
epipolaire 51 definie ci-avant. Plus precisement, autour de chacun de 
ces pixels selectionnes dans la zone de recherche, on genere (etape 
403) une deuxieme fenetre 61 g de meme caracteristique 
dimensionnelle que la fenetre cible. On stocke les niveaux de gris des 
pixels de chaque deuxieme fenetre et on procede (etape 404) k un 
traiternent de correlation entrc les niveaux de gris des pixels de chaque 
deuxieme fenetre 61g et les niveaux de gris des pixels de la fenetre 
cible 6 Id. 

Le traiternent de correlation permet done d'effectuer une mesure 
de similarite entre la fenetre cible et chacune des deuxiemes fenetres 
de la deuxieme image de fa9on k obtenir une valeur de similarite 
(valeur de correlation entre ces deux fenetres). 

C'est k partir de l'analyse de ces valeurs de correlation que Ton va 
proc6der k l'appariement de la region d'interet cible avec sa r6gion 
d'interet homologue. 

En ce qui concerne le traiternent de correlation (etape 404), on 
peut utiliser une correlation norm6e. Plus precisement, chaque valeur 
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de correlation d'une deuxieme fenetre sera alors donnee par la 
formule : 



ou N est le nombre de pixels d f une fenStre, k l'indice courant d'un 
pixel, I(k) le niveau de gris du pixel k de la fenetre cible et J(k)_le 
niveau de gris du pixel k de la deuxieme fenetre. 

On peut egalement utiliser une difference normee permettant 
d'obtenir la valeur de correlation par la formule : 



En theorie, la region d'interet homologue sera celle contenue dans 
la deuxieme fenetre qui presente la valeur de correlation maximale. 
Ceci etant, en pratique, plusieurs maxima de correlation peuvent etre 
trouves avec une intensite elevee. On selectionne alors (etape 405), 
parmi Tensemble des valeurs de correlation, un ensemble restreint de 
maxima de correlation selon un critere predetermine comme Tintensite 
du maximum ou comme le nombre de maxima dans la liste des valeurs 
de correlation rangees par ordre d'intensite decroissante. 

On peut egalement utiliser comme critere supplementaire de 
selection une determination de la dynamique des maxima selectionnes 
et la comparaison des valeurs de dynamique obtenues a un seuil 
predetermine. 

La determination de la dynamique de maxima est parfaitement 
connue de Thomme du metier. Cependant, celui-ci pourra 
eventuellement se reporter k la demande de brevet frangais n° 
91 15308 pour plus de details concernant une telle determination de 
dynamique. 

II convient de noter ici que le procede selon Tinvention permet de 
determiner et d'afficher sur Tecran du microordinateur, par exemple 



N 

^ I(k)J(k) 



I 
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par un clignotement, la region d'int^ret consideree comme etant la 
region homologue de la r6gion cible, et ce avec le minimum de risque 
d'erreur. ~._ " . . 

Ceci etant, en pratique, le procede selon Tinvention pourra 
selectionner un nombre restreint de regions d'interets candidates 
susceptibles d'etre homologues de la r6gion cible. Le radiologue aura 
alors la possibility de selectionner manuellement celle qu'il considere 
comme etant effectivement la region homologue. 

Ainsi qu'on Fa evoque ci-avant, on peut utiliser dans le traitement 
d'appariement, un traitement de comparaison specifique au lieu d'un 
traitement de correlation. 

Ainsi, comme illustre plus precisement sur la figure 4, on 
determine (etape 409) pour la region cible contenue a l'interieur de la 
fenetre cible 6 Id, et a partir des niveaux de gris des pixels contenus 
dans cette fenetre cible, une liste de caracteristiques permettant de 
caracteriser en general une region d'interet. Ces caracteristiques 
peuvent comporter par exemple une valeur representative de 
l'elongation de la region d'interet, une valeur representative de sa 
largeur moyenne, une largeur representative de son niveau de gris 
moyen, ou de son gradient moyen d'intensite. 

On effectue la meme operation (etape 410) pour les regions 
d'interets candidates contenues dans des deuxiemes fenetres 61 g 
definies a partir de pixels selectionnes dans la zone de recherche. A 
cet egard, afin de limiter le temps de traitement, les pixels 
selectionnes pour la determination des deuxiemes fenetres englobant 
les regions candidates, sont choisis parmi les pixels de la zone de 
recherche presentant des maxima de niveaux de gris dont la dynamique 
est superieure a un seuil predetermine. 

A partir de ces differentes listes de caracteristiques ainsi 
obtenues, on effectue un traitement de comparaison 411 entre la liste 
de caracteristiques affectees k la region cible et chacune des listes de 
caracteristiques affectees a chacune des regions candidates. 

En pratique, bien qu'on puisse envisager un autre type de 
traitement de comparaison, on utilisera de preference un traitement de 
comparaison prevoyant la minimisation d'une distance, par exemple la 
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distance eudidienne, entre la liste des caracteristiques de la region 
cible, consideree alors comme un vecteur, et une liste de 
caracteristiques d'une region candidate, consideree egalement comme 
un vecteur. 

On pourrait egalement prevoir, h. titre d'exemple, de determiner la 
valeur moyenne des valeurs absolues des differences entre les valeurs 
homologues de deux listes de caracteristiques, et de rechercher celle 
des valeurs moyennes qui est la plus faible pour en deduire la region 
d'interet homologue. 

La selection de la region homologue prevue dans Tetape 405 
s'obtient done par cette variante en minimisant une distance entre des 
listes de caracteristiques plutot qu'en maximisant une correlation 
normee ou une difference normee. 

L'invention prevoit encore certaines ameliorations. 

Ainsi, dans la variante de l'appariement prevoyant un traitement 
de correlation, on peut limiter la zone de recherche aux pixels de la 
zone epipolaire presentant des maxima de niveaux de gris dont la 
dynamique est superieure a un seuil predetermine. 

En outre, il est preferable, notamment pour s'affranchir des 
problemes de fond d'image, d'effectuer un filtrage prealable des 
images stereotaxiques du type dit "chapeau haut de forme" (ou 
transformation top-hat en langue anglaise). Ce genre de filtrage est 
connu de l'homme du metier et celui-ci pourra se referer pour plus de 
details k l'article de J. SERRA Image Analysis and Mathematical 
Morphology, Vol. 2, Academic Press 1988. 

Ce genre de filtrage permet d'extraire d'une image numerique les 
zones claires et etroites et permet de s'affranchir du fond. 

Les etapes du procede selon Tinvention qui vont maintenant etre 
decrites, et prevoyant notamment une reprojection sur Timage de 
centrage, s'apparentent ainsi notamment comme etant une verification 
des resultats de Tappariement precedent entre Timage cible et la 
deuxieme image, et perrnettent ainsi de minimiser encore les risques 
d'erreurs d'appariement. 

Ceci etant, bien que la suite du proc6de qui va maintenant etre 
decrite en se referant plus particuli&rement aux figures 4 et 6, prevoit 
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la possibilite d'utiliser des traitements de correlation ou de 
comparaison du type de ceux qui viennent d'etre d6crits, ce procede 
selon Tinvention utilisant une reprojection sur l'image de centrage, est 
independant de rappariement pr6alable utilise entre l'image cible et la 
deuxieme image. 

On suppose maintenant (figure 6) que le procede selon Tinvention 
a permis de sfilectionner sur l'image de gauche 12 deux regions 
d'interets candidates 5g et 20g susceptibles d'etre la region hornologue 
de la region cible 5d apparaissant sur l'image cible 10. 

A partir du pixel cible de la region d'interet cible et du pixel 
candidat ayant donne lieu k la selection de la region d'interet 
candidate 5g, on determine (etape 406), par un calcul trigonometrique 
de stereotaxic, et connaissant la g6ometrie du dispositif de stereotaxic, 
les coordonnees tridimensionnelles du point correspondant de l'objet 
d'interet candidat correspondant k ces deux regions 5d et 5g. Puis, par 
un calcul trigonometrique de stereotaxic inverse, on determine les 
coordonnees dans l'image de centrage du pixel projete correspondant k 
la projection du pixel de l'element candidat dans cette image de 
centrage. 

On remarque done que, dans le cas present, puisque la region 
candidate 5g est effectivement la region hornologue de la region 5d, la 
position du pixel projete dans l'image de centrage 1 1 correspond 
effectivement a une region d'interet projete 5p qui est effectivement la 
projection dans cette image de centrage de relernent d'interet 5, 

Par contre, par ces merries calculs, la region d'interet candidate 
20g va conduire k l'obtention d'une position projetee 20p dans l'image 
de centrage qui, dans le cas present, ne se superposera k aucune region 
d'interet physiquement presente dans cette image de centrage. 

Aussi, un traitement de correlation entre la fenetre cible contenant 
la region cible 5d et un voisinage du pixel projet6 contenant la region 
d'interet projetee 5p va donner une valeur de corr61ation sup^rieure k 
la valeur de correlation obtenue entre la fenetre cible 5d et un 
voisinage du pixel projete 20p. II en aurait 6t6 de mSme, meme si le 
voisinage du point projete 20p avait correspondu k une micro- 
calcification projet6e autre que celle correspondant effectivement k la 
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rnicrocalcification 5. 

Ce procede permet done de selectionner parmi les deux candidats 
5g et 20g la region 5g qui est effectivement la region homologue de la 
region 5d. 

Bien entendu, la conclusion aurait et6 la mSme en utilisant un 
traitement de comparaison 412 au lieu d f un traitement de correlation 
407. 

L'utilisation d'un voisinage du pixel projete, typiquement 
quelques pixels autour du point projete, permet de prendre en compte 
les erreurs eventuelles dues a un bouge de l'objet entre 1'exposition a 
0°et les expositions a ±15°. 

D'une fa9on plus generate, on peut ordonner l'ensemble des 
regions candidates obtenues dans la deuxieme image dans l'ordre 
decroissant de leur probabilite h etre la region homologue de la region 
cible. Ainsi, la region ayant par exemple donne le maximum de 
correlation avec la dynamique la plus forte sera affectee du rang 1 et 
ainsi de suite par ordre decroissant. On peut effectuer ce meme 
ordonnancement pour l'ensemble des regions projetees obtenues sur 
l'image de centrage. On selectionnera alors la region homologue de la 
region cible comme etant celle dont le produit des rangs respectifs 
dans les deux ordonnancements est minimum. 
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REVENDICATIONS 

1. Proced6 de localisation d'un element d'interet contenu dans un 
objet tridimensionnel a partir des positions de regions d'interets 
homologues correspondant audit element d'interSt et apparaissant dans 
un ensemble d'images st6r6otaxiques dudit objet, comportant une 
selection dans une premifere image sterdotaxique (10) d'une premiere 
region d'interet (5d) dite "cible", ainsi qu'un appariement de ladite 
premiere region avec une deuxieme region d'interSt (5g) homologue de 
la premiere et apparaissant dans une deuxieme image stereotaxique 
(12), caracterise par le fait que, les images stereotaxiques etant 
numerisees, on selectionne (400) dans Tetape de selection, un pixel 
cible (60d) dans ladite- region d'interet cible (5d), et dans l'etape 
d'appariement, on genere (401), autour du pixel cible selectionne, une 
fenetre cible (6 Id) de caract6ristiques dimensionnelles choisies et 
contenant ladite region d'interet cible, on determine (402) dans la 
deuxieme image, selon un critere de selection predetermine, un 
ensemble de pixels, on genere (403) autour de chaque pixel 
selectionne une deuxifeme fendtre de mSmes caracteristiques 
dimensionnelles que ladite fenetre cible, on effectue un traitement de 
correlation (404) entre les niveaux de gris des pixels de chaque 
deuxieme fenetre et les niveaux de gris des pixels de la fenetre cible, 
de fagon a obtenir une valeur de correlation pour chaque deuxieme 
fenetre, et on identifie la region d'interet homologue de la region 
d'interet cible a partir de I'analysc de Tensemble des valeurs de 
correlation ainsi obtenues, de fa?on k minimiser les risques d'erreur 
d'appariement entre les regions d'interet homologues. 

2. Procede selon la revendication 1, caracterise par le fait que 
l'analyse des valeurs de correlation obtenues comporte la selection 
d'un nombre choisi de maxima ou de minima de correlation, la 
selection de la region d'interet homologue s'effectuant parmi celles 
dont la valeur de correlation associ6e est l'un de ces maxima ou 
minima de correlation. 

3. Procede selon la revendication 2, caracterise par le fait que 
Tanalyse des valeurs de correlation obtenues comporte en outre une 
determination de la dynamique des maxima ou minima s61ectionn6s, et 
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la comparaison des valeurs de dynamique obtenues h un seuil. 

4. Procede selon Tune des revendications 1 & 3, caracterise par le 
fait que le traitement de correlation comporte un traitement de 
correlation normee. 

5. Procede selon Tune des revendications 1 & 3, caracteris6 par le 
fait que le traitement de correlation comporte un traitement de 
difference normee. 

6. Procede de localisation d'un el6ment d'interet contenu dans un 
objet tridimensionnel a partir des positions de regions d'interets 
homologues-correspondant audit element d'interet et apparaissant dans 
un ensemble d'images stereotaxiques dudit objet, comportant une 
selection dans une premiere image stereotaxique (10) d'une premiere 
region d'interet (5d) dite "cible", ainsi qu'un appariement de ladite 
premiere region avec une deuxieme r6gion d'interet (5g) homologue de 
la premiere et apparaissant dans une deuxieme image stereotaxique 
(12), caracterise par le fait que, les images stereotaxiques etant 
numerisees, on selectionne(400) dans l'etape de selection, un pixel 
cible dans ladite region d'interet cible, et dans l'etape d'appariement, 
on genere (401), autour du pixel cible selectionne, une fenetre cible 
(6 Id) de caracteristiques dimensionnelles choisies et contenant ladite 
region d'interet cible, on determine (409), pour ladite region d'interet 
cible, a partir des valeurs de niveaux de gris des pixels de ladite 
fenetre cible, un premier ensemble de valeurs numeriques 
respectivement representatives de caracteristiques predeterminees 
caracterisant toute region d'interet, et on stocke ce premier ensemble 
de valeurs numeriques, on determine (402) dans la deuxieme image, 
selon un critere de selection predetermine, un ensemble de pixels 
appartenant respectivement a des regions d'interet dites "candidates" 
apparaissant dans la deuxieme image, et on genere (403) autour de 
chacun de ces pixels selectionnes une deuxieme fenetre contenant 
ladite region d'interet candidate correspondante, on determine (410), 
pour chaque region d'interet candidate, a partir des valeurs de niveaux 
de gris des pixels de la deuxieme fenetre associee, un deuxieme 
ensemble de valeurs numeriques respectivement representatives 
desdites caracteristiques predeterminees, on stocke ce deuxieme 
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ensemble de valeurs numeriques, et on identifie la region d'interet 
homologue de la region d'interet cible k partir d'un traitement de 
comparaison (411) entre ledit premier ensemble de valeurs numeriques 
et chaque deuxifeme ensemble de valeurs numeriques, de fagon k 
minimiser les risques d'erreur d'appariement entre les regions d'interet 
homologues. 

7. Procede selon la revendication 6, caracterise par le fait que 
lesdites caracteristiques predetermines comprennent des 
caracteristiques de forme, de contraste, de gradient. 

8. Procede selon la revendication 6 ou 7, caracterise par le fait 
que le traitement de comparaison comporte un traitement de 
minimisation de distance entre les deux ensembles de valeurs 
numeriques. 

9. Procede selon Tune des revendications precedentes, caracterise 
par le fait qu'on determine dans la deuxieme image stereotaxique une 
zone epipolaire contenant au moins le segment epipolaire relatif au 
pixel cible, et par le fait que les pixels de la zone epipolaire 
contiennent les pixels selectionnes. 

10. Procede selon Tune des revendications precedentes, 
caracterise par le fait que lesdits pixels selectionnes sont ghoisis 
parmi les pixels de la deuxieme image presentant des maxima ou 
minima de niveaux de gris dont la dynamique est superieure k un seuil 
predetermine. 

11. Procede selon Tune des revendications precedentes, 
caracterise par le fait que l'appariement comporte prealablement un 
filtrage des images stereotaxiques, en particulier du type dit "chapeau 
haut de forme". 

12. Procede de localisation d'un 61ement dlnteret contenu dans un 
objet tridimensionnel k partir des positions de regions d'interets 
homologues correspondant audit 61ement d'interet et apparaissant dans 
un ensemble d'images stereotaxiques dudit objet, caracterise par le fait 
que, les images stereotaxiques 6tant numerisees, il comprend une 
selection dans une premiere image stereotaxique (10) d*une premiere 
r6gion d'interet (5d) dite "cible", une selection (405), dans une 
deuxieme image stereotaxique (12), k partir d'un premier appariernent 
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automatique, d'au moins une deuxieme region d'interet (5g) suceptible 
d'etre homologue de la region cible, la determination de la position 
spatiale d'un pixel candidat d'un element d'interet candidat 
correspondant k ces deux regions d'interet, la determination (406) de 
la position dans une troisi&me image stereotaxique, d'un pixel projete 
correspondant a la projection dans cette troisieme image stereotaxique 
du pixel candidat, et un deuxieme appariement automatique entre la 
region d'interet cible et un voisinage du pixel projete definissant une 
region d'interet projetee, de fagon a minimiser les risques d'erreur 
d'appariement entre les regions d'interet homologues. 

13. Procede selon la revendication 12, caracterise par le fait qu'on 
selectionne dans la deuxieme image stereotaxique, a partir dudit 
premier appariement automatique, un ensemble de deuxiemes regions 
susceptibles d'etre homologues de la region cible, on determine la 
position spatiale de chaque pixel candidat correspondant a chaque 
paire de regions d'interet formee par la region cible et Tune des 
deuxiemes regions, ainsi que la position dans la troisieme image de 
chaque pixel projete correspondant, et on effectue le deuxieme 
appariement automatique entre la region cible et un voisinage de 
chaque pixel projete. 

14. Procede selon la revendication 13, caracterise par le fait qu'on 
ordonne Tensemble des deuxiemes regions obtenues dans l'ordre 
decroissant de leur probability a etre la region homologue de la region 
cible, on ordonne I'ensemble des regions projetees obtenues dans 
l'ordre decroissant de leur probability a etre la region homologue de la 
region cible, et on selectionne la region homologue de la region cible 
comme etant celle dont le produit des rangs respectifs dans les deux 
ordonnancements est minimum. 

15. Procede selon Tune des revendications 12 a 14, caracterise par 
le fait que le premier appariement automatique est un appariement tel 
que defini dans Tune des revendications 1 a 11. 

16. Procede selon Tune des revendications 12 k 15, caracterise par 
le fait que le deuxieme appariement automatique comporte un 
traitement de correlation tel que defini dans Tune des revendications 1 
a 5, entre le voisinage de chaque pixel projete et une fenStre cible 
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contenant la region d'interet cible, ou un traitement de comparaison, 
tel que d6fini dans Tune des revendications 6 k 11 entre la region 
d'interSt cible et chaque region d'interet projetee. 
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